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O . ABSTRACT 

' Galactic stellar-population-synthcsis models, chemical-enrichment models, and long- 

duration Bulge microlensing events indicate about N tot = 10 s — 10 9 stellar-mass black 
CLi holes reside in our Galaxy. We study X-ray emission from accretion from the interstel- 

D ' lar medium on to isolated black holes. Although isolated black holes may be fewer in 

number than neutron stars, N^s ~ 10 9 , their higher masses, (M) ~ 9M Q , and smaller 
OO ' space velocities, a v ~ 40 km s _1 result in Bondi-Hoyle accretion rates ~ 4 x 10 3 times 

CnJ \ higher than for neutron stars. We estimate that ~ 10 4 isolated black holes within the 

Milky Way should accrete at M > 10 15 g s" 1 , comparable to accretion rates inferred 
for black- hole X-ray binaries given a total number of black holes N to t = iVglO 9 . If black 
*T . holes accrete at efficiencies only ~ 10~~ 4 (N^s /Ntot) ' 8 of the neutron-star accretion 

efficiency, a comparable number of each may be detectable. We make predictions for 
the number of isolated accreting black holes in our Galaxy which can be detected with 
. X-ray surveys as a function of efficiency, concluding that all-sky surveys at a depth 

of F = F_i 5 l(r 15 erg cm" 2 s" 1 dex" 1 can find N(> p) ~ 10 4 iV9(F_i 5 /e_ 5 )- L2 
isolated accreting black holes for a velocity dispersion of 40 km s _1 and X-ray accre- 
tion efficiency of e = e_5l0 -5 . Deeper surveys of the Galactic plane with Chandra or 
XMM-Newton may find tens of these objects per year, depending on the efficiency. We 
argue that a minimum mass can be derived for microlensing black-hole candidates if 
they are detected in the X-ray. 

O ■ 

!_i ' Key words: accretion, accretion discs — black hole physics — Galaxy: stellar content 

"o^ \ — X-rays: ISM, stars 
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1 INTRODUCTION 



5_i ' 

■ A black hole may only be detected via its interactions with surrounding matter or light. Black holes in binary-star systems are 
strong and transient X-ray sources fed by accretion from a stellar companion. However, isolated black holes — that is, black 
holes without companions or black holes in wide binaries — can easily escape detection. These may constitute the majority of 
the population of black holes since the fraction in close binaries with a < 100 AU is ~ 0.01 based stellar evolution calculations 
for massive star binaries (Chris Fryer, priv. comm. 2001). Young isolated neutron stars are much easier to detect as pulsars, 
but once the magnetic field, spin, and thermal energy decay away, an old neutron star becomes as invisible as a black hole. 
Several authors (Ostriker, Rees & Silk 1970, Treves & Colpi 1991, Blaes & Madau 1993) have proposed searching for isolated 
old neutron stars lit up by accreting matter from the interstellar medium (ISM). A tenth of the rest-mass energy of the 
accreting gas can be directly converted to luminosity assuming that all the accreted material is thermalised as it hits the 
neutron-star surface. These predictions motivated searches that turned up several isolated neutron-star candidates (Stocke et 
al. 1995, Walter et al. 1996, Wang 1997, Schwope et al. 1999); however, it is possible that these are young cooling neutron 
stars powered by processes other than accretion (Treves et al. 2000). 

The known Galactic black-hole X-ray-binary population has a velocity dispersion smaller by a factor of ~ 5 than neutron 
stars (Hansen & Phinney 1997, White & van Paradijs 1996, van Paradijs & White 1995) and average mass larger by a factor 
of ~ 6 (see below and Bailyn et al. 1998, Thorsett & Chakrabarty 1999). The accretion rate for an object of mass M accreting 
from the interstellar medium is determined by the Bondi-Hoyle formula (Bondi & Hoyle 1944), 

(«2 +c 2)3/2 ' W 
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where v is the velocity of the black hole with respect to the local interstellar medium, n is the number density of hydrogen 
in the interstellar medium, p, — p/n where p is the particle mass density, c s is the sound speed of the interstellar medium, 
and A is a parameter of order unity (we set A = 1 hereafter). Assuming that isolated black holes have the same mass range 
and velocity dispersion as black holes in X-ray binaries, one would expect the typical accretion rate on to isolated black holes 
from the ISM to be larger by a factor of ~ 4 x 10 3 compared to neutron stars. The expected detection rate of isolated black 
holes accreting from the ISM may be reduced as the number of black holes in our Galaxy may be smaller by a factor of ~ 10 
than the number of neutron stars, and the lack of a hard surface may result in lower luminosities for black holes accreting at 
small accretion rates. Garcia et al. (2001) argue that black holes have quiescent X-ray efficiencies ~ 10 -2 of neutron stars; 
however, Bildsten & Rutledge (2000) argue that the detected X-rays may have nothing to do with accretion, and V404 Cyg is 
an exception to this rule indicating that there may be a wide range of accretion efficiencies. If black holes do have a reduced 
efficiency compared to neutron stars, this reduction can be compensated by the larger accretion rates and the smaller scale 
height of black holes which increases their number density in the mid-plane where most of the interstellar gas mass is located. 
Black holes do not have an intrinsic magnetic field, so accretion can proceed uninhibited by magnetic pressure or torque. 
Thus, we are motivated to compute the expected numbers and accretion rates of isolated black holes in our Galaxy. 

The problem of black holes accreting from the interstellar medium has been considered by several authors. We approach 
the problem in a different way from previous authors (listed below) by using the properties of black-hole X-ray binaries and 
black-hole microlensing candidates to constrain the phase-space distribution of black holes within the entire Galactic disc 
(discussed in detail in Agol & Kamionkowski 2001, in preparation), while considering the properties of all phases of the ISM, 
deriving the distribution of accretion rates including radiative pre-heating, and finally deriving the distribution of X-ray fluxes. 
Grindlay (1978) considered spherical-accretion models for X-ray sources now known to be accreting from companions. Carr 
(1979) estimated the luminosity of black holes in the galactic disc, but did not predict their detection probability. McDowell 
(1985) computed the visual and infrared fluxes of black holes accreting spherically from molecular clouds, but only out to 
150 pc. Campana & Pardi (1993) computed the fluxes of black holes within the disc of the Galaxy accreting spherically from 
molecular clouds, without computing the probability distribution of accretion rates. Heckler & Kolb (1996) computed the 
fluxes of a putative halo population of black holes accreting spherically while passing within 10 pc of the Sun. The large 
velocities, low densities, and large scale heights of a halo population greatly reduces the number of black holes at a given 
accretion rate. Popov & Prokhorov (1998) took into account the spatial distribution of the warm interstellar medium, assuming 
the black holes receive a kick velocity at birth of ~ 200 km s _1 , which greatly reduces the accretion rates. Fujita et al. (1998) 
computed the fluxes of black holes in the disc accreting via an advection-dominated accretion flow at a distance of 400 pc (the 
Orion cloud) for a few velocities, without computing the full probability distribution. Grindlay et al. (2001) estimate that of 
order 10 3 black holes in molecular clouds might be detectable by the proposed Energetic X-ray Imaging Survey Telescope 
(EXIST) experiment based on an accretion efficiency of 10 -2 . 

The spectrum of an accreting black hole has been computed in the spherical-accretion limit (Ipser & Price 1977, 1982, 
1983); however, as we show below, the accreted material will possess non-zero angular momentum, possibly forming an 
accretion disc which may increase the efficiency of accretion. Thus, we parameterize our results in terms of the unknown 
accretion efficiency. 

In section 2 we summarize the known properties of binary and isolated black holes and their distribution in phase space. 
In section 3 we estimate the accreted angular momentum, discuss the range of accretion efficiencies, and discuss the effects 
of radiative feedback. In section 4 we discuss the properties of the interstellar medium necessary for computing the accretion 
rates of black holes. In section 5 we compute the distribution of accretion rates and luminosities for black holes within the 
disc of the galaxy. In section 6 we estimate the detection rates of black holes under the persistent and transient assumptions, 
and then we predict the fluxes as a function of distance for the three candidate black-hole microlenses. Finally in section 7 
we summarize. 



2 ISOLATED BLACK HOLES IN OUR GALAXY 

The relevant results from Agol & Kamionkowski (2001) on the number and distribution of black holes within the Galaxy are 
summarized here. 



2.1 Number of black holes 

Three microlensing events toward the Galactic Bulge have been found by the MACHO and OGLE teams (Bennett et al. 2001, 
Mao et al. 2001). These events are of long enough duration (> 1 year) to observe variations in the microlensing light curve 
due to the motion of the Earth around the Sun, an effect known as microlensing parallax. This breaks some of the degeneracy 
between distance, velocity, and mass of the lensing object; when combined with a model for the velocity and number density 
of black holes in the Galaxy, this results in an estimate of the likely mass of the lensing objects. In all three cases, the objects 
likely have masses greater than 3Mq, while they are not detected as main-sequence stars, leading to the conclusion that they 
are most likely black holes (Bennett et al. 2001, Agol & Kamionkowski 2001). Based on the probability of microlensing toward 
the Bulge and assuming that isolated black holes have the same scale height as black holes in binaries, one concludes that 
the total number of black holes in the Milky Way disk is 2 x 10 8 < 7V to t(//0.2) < 4 x 10 9 at 95 per cent confidence, for a 
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microlensing detection efficiency /, estimated to be 20 per cent (Alcock et al. 2000). Given the uncertainty in these estimates, 
we will make predictions based on TVtot = A9IO 9 black holes, and note that our results scale with Ng. 

As discussed below, there is some evidence that the velocity dispersion of isolated black holes is similar to that of the 
stellar disk. We thus assume that the number density of black holes is distributed as 



Nbh = N® exp 



\A 

H 



(R-R Q 



Rdis 



(2) 



where Nq is the number density at the solar circle, Z and R are the cylindrical Galactic coordinates, Rq = 8 kpc is the 
distance to Galactic centre, H and -Rdisc are the disc scale height and scale length. We will assume that the black-hole disc scale 
height is H = 375 pc and scale length is -Rdisc = 3 kpc. This implies a local number density of Nq = 1.64x 10 5 -/Ve,(ii/y375pc) _1 . 



2.2 Velocity distribution of isolated black holes 

White & van Paradijs (1996) have estimated that the rms distance from the Galactic plane for black-hole X-ray binaries, 
z r ms, is about 450 pc, which corresponds to a scale height of H ~ 320 pc and velocity dispersion of a v ~ 40 km s . This is 
consistent with the radial-velocity measurements of low-mass black-hole X-ray binaries when corrected for Galactic rotation 
(Brandt et al. 1995, Nelemans et al. 1999), with the exception of GRS 1655-40 which has v r — —114 km s _1 . We have 
improved on the work by White & van Paradijs (1996) by using a maximum-likelihood analysis based on 19 black-hole X-ray 
binaries, concluding that H = 375lgQ 5 P c (Agol & Kamionkowski 2001). In addition, the three isolated black holes have a 
vertical velocity dispersion of a z < 27 km s _1 at 95 per cent confidence, indicating that their kinematics, and thus spatial 
distribution, is similar to that of black holes in X-ray binaries. Thus, we take the isolated black holes to have the same velocity 
dispersion and scale height as black hole X-ray binaries. 



2.3 Mass distribution of isolated black holes 

The distribution of masses of black holes in X-ray binaries and microlensing events can be fit with a power law mass distri- 
bution, 

dN/dM = n M M" 7 , (3) 

with a lower mass cutoff at Mi and upper mass cutoff of M2. We obtain a best fit for 7 = 0.6, Mi = 4Mq, and Mi = 14Mq 
(Agol & Kamionkowski 2001). Since there may be larger-mass black holes that have not yet been detected, we impose an 
upper limit mass cutoff of M2 = 5OM0 based on the theoretical mass limit for pulsational stability (Schwarzchild & Harm 
1959), assuming that about half of the mass is lost before collapse. This yields Mi = 4M© and 7 = 2.4, close to the value 
of the Salpeter IMF (Agol & Kamionkowski 2001). The average mass of either distribution is ~ 9M0, close to the mean of 
the data (M) = 9.6M©. We have normalized this to unity, um = (7 — l)/(A4l~ J — A/j -7 ), and we assume that it applies 
throughout the Galaxy. We will later mention how our predictions differ for M2 = 13M© , but will otherwise use M2 = 5OM0 . 



3 EFFICIENCY OF ACCRETION 

3.1 Angular momentum of accreted material 

The Bondi-Hoyle formula may be appropriate for computing the accretion rate, but not the efficiency. If the accreting gas 
were perfectly uniform, then the accreted angular momentum is quite small, and the accretion flow follows a one-dimensional 
supersonic solution, for which the luminosity can be predicted (Shapiro & Teukolsky 1983, Ipser & Price 1977). However, small 
perturbations in the density or velocity of the accreting gas will lead to an angular momentum large enough to circularise gas 
before it accretes (Shapiro & Lightman 1976). If there is a gradient in the density field that is perpendicular to the direction 
of motion of the black hole, then the accreted angular momentum scales as 

l=\^vr A , (4) 
4 p 

where Ap is the difference in velocity between the top and bottom of the accretion cylinder and ta = GM/(v 2 + (? s ) is the 
accretion radius. Numerical simulations confirm that this formula is correct to order of magnitude (Ruffert 1999). The observed 
density fluctuations in our galaxy scale as Sp/p ~ (L/10 18 cm) 1,/3 (Armstrong, Rickett, & Spangler 1995), extending down to 
a scale of ~ 10 8 cm. The exponent of this scaling is nearly consistent with a Kolmogorov spectrum of density (Dubinski et al. 
1995). Evaluating this at L = 2rA, we can then find the radius of the resulting accretion disc, rdisc, by equating the angular 
momentum of the gas with the Keplerian angular momentum due to the black hole, ?k cp = VGMrdisc- This gives 

™==10*r ( M\ 2/3 ( V^3 ^ 10/3 
r g 9 \9M Q J ^40kms- ] 



(•») 
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So, we see that a disc will almost always form in interstellar accretion. A similar argument is given in Fujita et al. (1998) for 
the case of molecular clouds, for which the slope of velocity fluctuations is similar to that expected by a Kolmogorov spectrum 
(Larson 1981). The neutral interstellar medium seems to follow a power-law with a similar slope (Lazarian & Pogosyan 2000), 
which will also lead to non-zero accreted angular momentum; however, the length scales probed in molecular and atomic gas 
are much larger than the scale of the accretion radius. 

The accretion time scale for the disk scales as td = a^g(GM/r| isc )~ 1 ' /2 (rdi S c/ft) 2 where h is the height of the disc and 
ass is the ratio of the viscous stress to the pressure in the disk. If the velocity and density fields of the ISM are incoherent on 
the scale of these fluctuations, then the angular momentum can change on a time scale of t a = ta/v as gas will be accreted 
with the angular momentum vector pointing in a different direction after this time scale. Thus, the angular momentum of the 
disc can be reduced due to a random walk of the accreted angular momentum vector if t a <$Ltd- We estimate 



which indicates that there may only be a slight reduction in the angular momentum since these time scales are comparable. 
The residual angular momentum means that we cannot use spherical-accretion formulae for the luminosity, but we must 
estimate the efficiency of the accretion disc that forms. 

3.2 Range of efficiencies 

We assume that the black holes radiate a spectrum with spectral index a = 1, where F v oc v~ a , consistent with the typical 
spectrum of black-hole X-ray binaries, so that equal energy is emitted per unit decade of energy. We define L[ on as the 
luminosity from 1 to 1000 Ryd (13.6 eV), so the efficiency is ei on = £ion/(Mc 2 ). Estimates of the accretion rate in Cygnus X-l 
indicate that it is accreting with an X-ray efficiency of £i on ~ 10 _1 (Shapiro & Teukolsky 1983). However, at low accretion 
rates the density of the gas decreases, so it may not be able to cool, reducing the radiative efficiency. By estimating the mass 
accretion rates in black-hole X-ray binaries and X-ray luminosities in the quiescent state, Lasota (2000) concludes that the 
2-10 keV efficiency is about 2 x 10 -6 assuming that this luminosity is solely due to accretion, implying eion ~ 10 -5 . Low 
upper limits on the X-ray Bondi-accretion efficiencies are inferred for some supermassive black holes in elliptical galaxies 
(Loewenstein et al. 2001). The Galactic-centre black hole has an X- ray efficiency of ~ 4 x 10 -8 - 2 x 10~ 6 in the 2-10 keV 
band for a Bondi-Hoyle accretion rate of 1O _6 M0 yr _1 (Baganoff et al. 2001). Thus, we parameterize our results in terms of 
an unknown X-ray efficiency. 

For certain ranges of parameters, the disc may be unstable to the hydrogen- ionisation disc instability (Lasota 2001), 
creating intermittent outbursts of higher luminosity between long periods of quiescence, with changes in luminosity of 10 5 . 
We will discuss the possible effect this may have on the number of detectable sources in section 6.2. 

3.3 Radiative feedback 

The gas in the vicinity of the black hole will be heated by the ionising radiation field produced by accretion. If the temperature 
rises enough, then the accretion radius, and thus the accretion rate, will decrease (Shvartsman 1971, Ostriker et al. 1976). For 
an a = 1 spectrum, the mean opacity from 1-1000 Ryd of neutral gas with cosmic abundance is <7 C oid = 3.87 x 10 -18 cm 2 per 
hydrogen atom. The heating rate is then given by H = L x nHO~ CO M/ '(47r_R 2 ), where R is the distance from the black hole. For 
temperatures smaller than ~ 10 4 K and for a small ratio of radiation energy density to gas energy density, the heating rate 
greatly exceeds the cooling rate. Consider a fluid element directly in front of the path of a moving black hole. As the black 
hole approaches with velocity v, the ratio of the total heating rate to the initial energy density, eo of the gas is 



where e to t is the total energy density absorbed by a fluid element (neglecting cooling), v = «4o40km s _1 and T — T4IO 4 K. 
Thus, even for small efficiencies the gas is strongly heated outside the accretion radius. Once H and He are ionized, the mean 
opacity drops by three orders of magnitude, and the gas remains at around 2 x 10 4 K. This has been confirmed with detailed 
numerical simulations by Blaes, Warren & Madau (1995). Above this temperature, the cooling and heating times become 
much longer than the accretion time so the gas falls in adiabatically. Thus we take the temperature of the accreted gas to be 
To = 2 x 10 4 K for all phases of the ISM except the hot HII. 

4 GAS IN OUR GALAXY 
4.1 Phases of the ISM 

The interstellar medium in our Galaxy consists of at least five identifiable phases (Bland-Hawthorn & Reynolds 2000) in 
approximate pressure equilibrium, P ~ 2 x 10 3 fcs. The highest accretion rates will occur in the regions of the highest density. 
This means that the observable sources will be dominated by the giant molecular clouds (GMCs) consisting mostly of H2 




(6) 




(7) 
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(T ~ 10 K) and the cold neutral medium consisting of HI clouds (CNM, T ~ 10 2 K). These phases comprise approximately 
40 and 50 per cent of the total mass of gas in our Galaxy. 

The warm neutral and ionised components (T ~ 8000K) of the interstellar medium have densities that are ~ 10" 2 of 
the CNM, so the maximum accretion rates are ~ 10 -2 times smaller than than the maximum for the CNM, while the hot 
component of the ISM (T ~ 10 6 K) has a density ~ 10~ 2 of the warm components, reducing the maximum accretion rate 
by another factor of 10 -2 . However, these phases fill a much larger volume in the interstellar medium, and thus should be 
considered. 

We next refine the above estimates by considering the distribution of each phase as a function of number density. 



4.2 Filling fraction of interstellar medium 

To estimate the distribution of accretion rates for black holes in our Galaxy requires knowing the volume of the Galaxy filled 
by gas of a given number density. We define this as 

£=/.»-», (8) 

where df/dn is the fraction of volume filled by gas with number density between n and n + dn, and fo is a normalization 
constant. Berkhuijsen (1999) has shown that (3gmc ~ 2.8, while Pcnm ~ 3.8. This can be seen for the molecular clouds by 
noting that the mass function of molecular clouds in our Galaxy scales as dN/dM oc M -1 ' 6 (Dickey & Garwood 1989), while 
number density of clouds scales with size, D, as n oc D^ 1 (Larson 1981, Scoville & Sanders 1987). Assuming that clouds have 
a spherical volume, one finds df/dn oc V (dM / dn) (dN/dM) oc n~ 28 . Since the number density-size relation is derived for the 
average density of clouds, /3 may be slightly smaller since each cloud may contain regions where the gas density is higher than 
average. 

We assume that the GMCs span number densities from ni = 10 2 to 112 = 10 5 cm~ 3 , while the CNM spans number 
densities from ni = 10 to ni — 10 2 cm" 3 . Finally, we assume that the scale height of the molecular clouds is 75 pc, while the 
scale height of the CNM is 150 pc. 

The surface mass density of molecular gas as a function of radius we take from Clemens, Sanders & Scoville (1988) and 
Scoville & Sanders (1987), while the distribution of the CNM we assume is constant with E = 4.5M© pc -2 for R > 4 kpc, 
and zero inside 4 kpc (Scoville & Sanders 1987). The molecular-gas surface density has a a strong peak at the Galactic centre, 
another peak around 5 kpc from the Galactic centre, then decreases outwards (Scoville & Sanders 1987). We ignore variations 
of the surface density in azimuth, and we approximate the distribution in z as an exponential, fo(z) = /o(0) exp(— \z\/H). 
From these surface densities, we compute the filling fraction in the mid-plane as, 



S(/3-2)(ra^-n^) 
2^-l)(n^-n^)' 



M0) - — - — — 2 _ 2^37, (9) 



where a = 2.72m p for the GMC, and a = 1.36m p for the CNM. At the solar circle, this gives /o(GMC) = 10~ 3 , while 
/o(CNM) = 0.04. 

The warm HI, warm HII, and hot HII we take as constant densities of 0.3, 0.15, and 0.002 cm -3 respectively, mid-plane 
filling factors of 35, 20, and 40 per cent, as well as scale heights of 0.5, 1, and 3 kpc (Bland-Hawthorn & Reynolds 2000). We 
assume that the hot HII has a temperature of 10 6 K, and thus c s = 150 km s _1 . The corresponding surface mass densities 
are 4, 2, and 0.2 Mq pc -2 , respectively, which we assume are constant as a function of Galactic radius. 



5 ACCRETION RATE AND LUMINOSITY DISTRIBUTION FUNCTIONS 

The distribution function of mass-accretion rate for black holes located at a given region in the Galaxy is given by 

M2 n 2 00 



j£N_ 

dMdV 



= dM dn 



Mi ni 

where dV is the volume element and 



^^=iVBH(i?, 2 )W-^exp 



V 



n M M~ 



(10) 



(11) 



We define M = 7rG 2 M 2 5 m p (lkm s" 1 )" 3 = 3.8 x 10 14 , so that M(n, M,v) = M nM 2 (v 2 + c 2 )~ 3/2 . For a given mass and sound 
speed, we can define the minimum number density required for accretion at a rate greater than M: n > no — Mc 3 s /(MqM 2 ). 
Using equations (|l|) and (^), we can first carry out the v integration analytically, 



M 2 



max no ,nn 



j£N_ 

dMdV 



N B 



dM n M M~ 



, . [2v (MonM 2 ) 2 ' 3 

dHfon MV3 6XP 



Ml 



max{ni ,tiq } 



2a 2 



(12) 
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where v 2 = (M nM 2 /M) 2/s - c 2 . The remaining two integrals we compute numerically. 
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Fig. 1: The number density of black holes accreting with accretion rates between M and M + dM as a function of M at 
the solar radius for a v = 40 km s _1 . The solid curve refers to black holes accreting within molecular clouds, while the short 
dashed curve refers to black holes within the cold neutral medium. The dotted curve is black holes within the warm HI, the 
dash-dot curve for the warm HII, while the dash-triple dot curve for the hot HII. 

In Figure 1 we plot the function dN/dM for the various phases of the interstellar medium at the solar circle. The densest 
gas (GMCs) dominates the highest accretion rates, while the hottest gas dominates the lowest accretion rates. For the hot 



HII, accretion is subsonic and n is assumed to be constant, so M oc M . Thus, dN/dM oc M 



"(l+7)/2 



spanning a range in M 



of (Mi/M\) , or two decades. This is consistent with the numerical spectrum shown in Figure 1. The other phases have more 
complicated accretion-rate distributions since a v > c s . 

To compute the expected total number of black holes in the Galaxy accreting at a rate M, we integrate the luminosity 
functions (assumed to be constant) over the gas filling fraction times the number density of black holes as a function of 
position over the entire Galaxy, 



d 2 N 



dN _ 
dM ~ dMdV 



dR 2nR exp 



R — Ri 



o 



Rd 



(fl o ,z=0) o 



J E(R e ) 



dz exp 



-HBH Hgas 



(13) 



Figure 2(a) shows N(> M) = J- dM' (dN/dM 1 ) for both black holes and neutron stars. For neutron stars, we choose 

700 km s" 1 (Cordes & Chernoff 1998), 



per cent with a v 

-3 



175 km s 



and 14 per cent with a v 

-,9 



M = 1.4M Q , 8 

iV© = 5.2 x 10 5 kpc -3 , and H = 1 kpc, yielding a total of 10 a neutron stars within our galaxy. The total number of accreting 
black holes is smaller than the total number in the disk since we have assumed that the gas filling factor decreases with scale 
height. This reduction may be artificial since gas at some density presumably fills the entire Galaxy, but this only affects 
small accretion rates. 

We expect ~ 10 4 iVg isolated black holes to be accreting at rates comparable to the estimated rates in black-hole X-ray 
binaries, M > 10 15 g s _1 (van Paradijs 1996). This exceeds the parent population of black- hole X-ray binaries, which may be 
10 2-3 (Iben et al. 1995). Black holes outnumber neutron stars due to their larger mass and smaller space velocity as shown 
in Figure 2(a). 

To compute the number of detectable isolated black holes requires converting the accretion-rate distribution into a 
luminosity function. Given the uncertainties in accretion physics, we will make two assumptions based on what is known 
about black-hole X-ray binaries: persistent and transient sources. The first assumption shown in Figure 2(b) is that the 
sources are time-steady, similar to high-mass X-ray binaries such as Cygnus X-l. Under this assumption, we simply convert 
the accretion rate into an X-ray luminosity in a given X-ray band. Shown is the ionising radiation band from 1-1000 Ryd for 
Eion = 10 -5 ; the horizontal axis simply scales with efficiency. The black-hole and neutron-star functions overlap for ej on ~ 10~ 5 
despite the lower black-hole efficiency due to the larger black-hole mass, smaller velocity dispersion, and smaller scale height. 
The total luminosity of isolated accreting black holes within the Milky Way is ~ 10 36 A'9e_5 erg s _1 . 
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Fig. 2: (a) The number of black holes in the Milky Way with accretion rates greater than M for a v = 40 km s 1 (solid line). 
For comparison, the dotted line shows the same computation for neutron stars, (b) The number of black holes in the Milky 
Way with luminosities greater than Lion with ei on = 10 -5 and the number of neutron stars for ei on = 10 -1 (dotted line). 

The second assumption that the black holes are transient we discuss in section 6.2 below. 



6 DETECTION PROBABILITY 

To compute the detection probability, we need to integrate the number of sources above a given observed flux as a function of 
position within the Galaxy. We will consider two energy emission bands, 'soft' X-ray from 1-10 keV, and 'hard' X-ray from 
10-100 keV, which have efficiencies of e so ft and e nar d. Since the model spectrum we have used is flat, uf v oc v , the luminosity 
per unit decade of photon energy is constant. Thus, the flux in each of these bands is 1/3 of the flux of ionising photons which 
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covers three decades, or e so ft = £hard = £ion/3. A more complex X-ray spectrum will change our results by factors of order 
unity, which we defer to future work given the poor theoretical understanding of the spectra of accreting black holes. 

For neutron stars, we assume a spectrum composed of a blackbody with T = 0.15 keV plus a power law with a — 1 
(i.e. flat in vF v ) and the flux from 1-10 keV equal to 1 per cent of the total flux of the blackbody component. This was 
chosen to be similar to the spectra of neutron-star binaries Cen X-4 and Aquila X-l in quiescence (Rutledge et al. 2001a,b). 
This spectrum yields e so ft = 0.13ei on and thard = 0.03ei on due to the fact that most of the flux is in the extreme ultraviolet. 
We neglect photoelectric absorption since the half-column density of molecular clouds is approximately Nh ~ 1 x 10 22 cm -2 , 
which leads to absorption of only 36 per cent of the flux in the soft X-ray band (using the opacities of Morrison & McCammon 
1987). Toward the Galactic centre, the hydrogen column density is about Nh ~ 6 x 10 22 cm -2 due to the ring of gas near the 
centre, leading to absorption of about 86 per cent of the soft flux. In the hard-X-ray band, the absorption is negligible. 



6.1 Persistent sources 

For the persistent-source assumption, we assume that the accretion flow emits isotropically, so we simply need to integrate 
the luminosity function over the entire galaxy with the transformation: 

dN dN AtvD 2 , . 
= (14) 

dFdV dMdV Esoft, hardC 2 ' 

where D is the distance to a given point in the Galaxy. Figure 3 shows the predicted number of sources as a function of flux 
integrated over the entire Galaxy for the persistent-source assumption. Despite the much lower assumed efficiency, the number 
of detectable black holes becomes comparable to the number of detectable neutron stars for €bh ~ W~ 4 £ns(Nbb /Nns) ' 8 ■ 
The number of black holes above a certain flux at high flux scales approximately as N(> F) — 10 4 F -1 ' 2 near F = 10~ 15 e_5 
erg cm -2 s _1 . Given the flatness of this dependence, the best detection strategy is to cover as much area of the sky as possible, 
assuming the constraint of a fixed amount of observing time with an X-ray telescope with a sensitivity that scales as t -1 ' 2 . 
Reducing the mass cutoff from — 50M Q to M2 = 13Mq only decreases the number of objects by a factor of ~ 2 since 
there are fewer objects at high mass (Figure 3). 




Fig. 3: The number of black holes visible to an observer located within the Solar system as a function of X-ray flux in either 
the soft or hard X-ray bands for a persistent accretion efficiency of ejon = 10 _J with e S oft,hard = £km/3 for M2 = 5OM0 (solid 
line) and M2 = 13Mq (dashed line). For comparison the number of neutron stars above a given soft X-ray flux is shown as 
well in the soft X-ray band with e S oft = 0.13ei on = 0.013 (dotted line). In addition, we show predictions for a population of 
10 6 black holes with mass 25OM0 for efficiencies of 10 -5 (dash-dot line). Note that the horizontal axis scales with efficiency 
for each object. 

We have also computed the the number density of black holes with fluxes greater than 10 _15 e_5 erg cm -2 s _1 as a function 
of Galactic latitude and longitude for one quadrant of the Galaxy (Figure 4) . The black holes are strongly concentrated toward 
the Galactic plane and toward the Galactic centre. Half of the black holes lie within the dashed contour which covers an area 
of 504 deg 2 in Figure 4, giving an average of 3Ag per square degree within this region near the Galactic plane. We have not 
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attempted to include the inhomogeneity of the ISM which might be important for nearby black holes in the local bubble or 
for black holes in interstellar clouds. 
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Fig. 4: The logarithm of the number density of black holes per degree square with fluxes above Fx = 10~ 15 e_5 erg cm -2 s" 1 
in the soft or hard band as a function of galactic longitude, and latitude, b(°). The labels indicate the logarithm of the 
number density. The dashed line indicates the contour within which 50 per cent of the black holes are contained. The total 
number of black holes above this flux level is ~ 5000, while the maximum number density is 94 per square degree at the 
Galactic Center. 

Finally, we have computed the distribution of sources with distance from the Sun. Figure 5 shows the number of black-hole 
sources with fluxes greater than F — 10~ 15 e_5 erg cm -2 s" 1 less than a given distance for each of the phases of the ISM. The 
total number of sources above this flux level is 38OOA9 in the GMCs, 67007V9 in the CNM, 900iVg in the warm HI, 200A r 9 in 
the warm HII, and in the hot HII. The average distance for the detectable black holes is 5 kpc in the GMCs, 1 kpc in the 
CNM, 200 pc in the warm HI, and 140 pc in the warm HII. The average distance of black holes in all phases is 2 kpc. Note 
that there is a rise in the number of black holes in GMCs at the Galactic centre (8 kpc away) since the surface density of 
molecular gas increases by a factor of ~ 10 within 500 pc of the Galactic centre. The sources in the gas clouds tend to be 
more distant since it requires looking further in the Galaxy to find clouds, while the the sources in the interstellar gas tend 
to be closer since they are fainter due to lower accretion rates. 

The sensitivity of the ROSAT All-Sky Survey (RASS) Bright Source Catalog is about 10" 12 erg cm" 2 s" 1 (Voges et al. 
1999; we have converted count rate to flux assuming a = 1), so at most 3Ngellc, accreting black hole should have been detected 
by this survey given the assumption of persistence. We can see immediately from Figure 3 why it may have been difficult to 
detect isolated accreting neutron stars with the RASS since we predict about one visible neutron star on the sky at this flux 
level, which will be reduced further by X-ray absorption. Danner (1998) has surveyed sources in molecular clouds within the 
RASS, finding no new isolated neutron-star candidates, and thus, we presume, no isolated black-hole candidates. For an X-ray 
observatory such as EXIST, which will reach sensitivities of ~ 10" 13 erg cm" 2 s" 1 for one-year co-added data (Grindlay et al. 
1999 ), about ~ SONgells sources may be detectable at hard X-ray energies as well. EXIST will have an angular resolution of 
about 2 arcminutes and position sensitivity of ~ 30 arcsec, which should make identification feasible within molecular clouds. 

The Chandra and XMM-Newton observatories have much better sensitivities ~ 10~ 14 erg cm" 2 s" 1 for a 10 kilo-second 
(ksec) observation; however, both satellites cover a much smaller field of view (16' x 16'). The ChaMPlane survey (Chandra 
Multi-wavelength Galactic Plane Survey, Wilkes et al. 2000) will cover about 1 square degree per year to a flux of ~ 10~ 15 
erg cm" 2 s" 1 for the 1-10 keV band. This flux limit is deep enough that ~ 10 4 A r geL| sources within the Galactic plane might 
be detectable (Figure 3); however, the central 50 per cent are spread over an area of 504 deg 2 in the Galactic Plane (Figure 
4). The XMM-Newton serendipitous Survey will cover an area of the sky ~ 10 times larger than the Chandra Serendipitous 
Survey at similar sensitivity (Watson et al. 2001). Thus, these surveys may result in ~ SONgel:^ detections per year. 

One might expect that targeting molecular clouds may be advantageous for finding black holes. However, due to their 
proximity, the 23 clouds within 1 kpc have angular sizes of about 5-20° (Dame et al. 1987), which will make them rather 
difficult to cover with the small field of view of Chandra or XMM-Newton. For instance, the single pointing in the p-Ophiucus 
molecular cloud observed by Imanishi et al. (2001) covers a volume ~ 18 pc 3 (assuming the cloud is spherical), so we expect 
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only ~ 3 x 1Q~ 2 Nq black holes to reside within the observed volume. Observations pointed at molecular clouds will not 
improve these detection rates, but observations pointing toward the Galactic centre should improve the number of detections 
significantly. Surveys targeted at the Galactic Center (Baganoff et al. 2001) may find 7N^e]l1 per Chandra or XMM-Newton 
field. 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

log(D(kpc)) 

Fig. 5: The total number of black hole sources with fluxes greater than F = 10 _1 e_s erg cm - s _1 and less distant than D. 
The solid line is for the CNM, the dotted line for the GMCs, the dashed line for the warm HI, and the dash-dot line for the 
warm HII. 



6.2 Transient sources 

The second assumption is that each of the sources are transient, driven by the hydrogen-ionisation instability at radii in the 
accretion disc with T ~ 5000 K as in low-mass X-ray binaries. In this case, we assume that the sources undergo outbursts 
for several months with luminosities near Eddington, while remaining faint for several decades. Chen, Shrader & Livio (1997) 
have shown that soft X-ray transients tend to peak at L out ~ 0.2Z/Edd with an average duration (decay time scale) of t ~ 20 
days and recurrence times varying from 2 years to 60 years for black holes. These are only the observed transients; transients 
with longer recurrence time scales may not have had time to repeat. If we assume that each black hole accretes quiescently 
for a time t q , storing up mass in the accretion disk with a small fraction of gas accreting on to the black hole, and then all of 
this mass is released in an outburst of duration t , then 
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where e is the radiative efficiency of accretion during outburst. The outburst rate in the Milky Way is then 
N out = / dM 

J Mi JO 



'j.v d2N - 
dM -t a 

dMdM q 



(15) 



(16) 



Assuming a 10 per cent efficiency during outburst, we find Nout ~ 536A9 yr~ for a v — 40 km s _1 . This is clearly inconsistent 
with observations, demonstrating that not all isolated black holes undergo transient outbursts if the assumptions of our 
calculation are correct. If only black holes accreting above M > 10 15 g s _1 undergo outbursts, consistent with the observed 
X-ray novae (van Paradijs 1996), then these numbers reduce to 50 yr _1 , still large given that about ~ 1 transient is detected 
each year, and generally they are found to have companions. We conclude that only a small fraction of isolated black holes 
might experience transience similar to X-ray novae. 



6.3 Detecting microlensing candidates 

The uncertainty in distance of the MACHO black-hole microlensing candidates leads to an uncertainty in their mass, and thus 
an uncertainty in their status as black holes. One possible way to constrain the mass of these lenses is to look for X-rays from 
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accretion. Figure 6 shows the flux as a function of distance for MACHO 96-BLG-5, MACHO 98-BLG-6, and OGLE-1999- 
BUL-32 (Mao et al. 2001) for a number density of 1 cm -3 with the total velocity set to times the observed sky velocity 

projected on to the lens plane (ignoring the unknown projected source velocity) and subtracted from the Galactic rotation 
velocity (220 km s" 1 ) and an X-ray efficiency of e\ on = 10~ 3 . We note that if the sky velocity of the source is significant, then 
the inferred lens velocity will be different, typically smaller, which will increase the luminosity. The flux scales as F oc enD~ A 
for small distances since the mass depends on the inverse of the distance and the flux decreases as the inverse of distance 
squared. Observing in the X-ray would be advantageous to avoid confusion with the background microlensed star. In the soft 
X-ray, a 1-10 keV flux of ~ 10 -15 erg cm -2 s _1 can be detected with a ~ 400 ksec exposure with Chandra or ~ 200 ksec with 
XMM-Newton. A detection of a black-hole candidate at such a flux level will allow an upper limit on the distance of ~ 1 — 3 
kpc to be placed, assuming they accrete from the diffuse ISM with n < n max = 1 cm" 3 and assuming an X-ray efficiency less 
than e max = 10~ 3 . A decrease in efficiency will lead to a decrease in flux, and thus a decrease in the distance upper limit. An 
upper limit on distance converts into a lower limit on mass, since M oc D~ , which may confirm the black-hole hypothesis. 
The mass lower limit is proportional to M m j n oc (tmaxnmai/f) -1 ' 4 , a weak dependence on our assumed parameters. A cloud 
or companion could increase the accretion rate, resulting in a higher luminosity, thus a larger distance for the same flux; 
however, clouds occupy less than 5 per cent of the volume near the Galactic plane and can be searched for in HI 21 cm or 
CO absorption or emission (the MACHO fields are toward Baade's window, thus unlikely to contain a high molecular column 
density), while companions should be detectable with HST. 



-10 




log(D(kpc)) 

Fig. 6: The flux as a function of distance for MACHO 96-BLG-5 (solid), MACHO 98-BLG-6 (dashed) and OGLE-1999-BUL-32 
(dotted) as a function of distance for ej on = 10 -3 . 



7 CONCLUSIONS 

We have estimated the number of isolated black holes that might be revealed by their X-ray emission from accretion of 
interstellar gas. We have improved upon previous calculations by taking into account the density distribution of the interstellar 
medium, as well as using the known properties of black-hole X-ray binaries and MACHO black-hole candidates to constrain 
the phase space and accretion efficiency of isolated accreting black holes. We conclude that persistent isolated black holes may 
be competitive with neutron stars in creating detectable X-ray flux if their efficiencies differ by cbh ~ W 4 £ns(Nbh /Nns) ' 8 ■ 
The ROSAT survey did not have the sensitivity to detect either isolated accreting neutron stars or black holes; an all-sky 

~\a 2 1 3 12 

survey with two orders of magnitude more sensitivity, ~ 10 erg cm s dex— 1, may be able to detect ~ 10 Nae^s 
accreting black holes. The X-ray observatories XMM-Newton and Chandra may be able to detect this population with pointed 
observations; however, they suffer from having small fields of view. Despite this limitation, Galactic plane surveys being carried 
out with these telescopes may be able to detect tens per year for e_5 = 1. If the accreting sources with M > 10 15 g s" 1 
go through transient outbursts, then 50 isolated black holes per year would be detectable as X-ray novae, so we reject this 
possibility since isolated X-ray novae have not been observed. 

The greatest uncertainty in our predictions is the efficiency of black-hole accretion, with which we have parameterized 
our results. Most of the detectable black holes should reside in interstellar clouds that have higher densities; however, this 
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leads to the problem of confusion with other X-ray sources such as the coronae of massive stars. Some tests for whether 
an accreting object is a black hole are the following: 1) Is there high-energy emission? Accreting black holes tend to show 
spectra which have power laws extending up to ~ 10 2 keV (Grove et al. 1998). 2) What is the nature of the variability? 
Accreting black holes show no pulsations, show only QPOs with v < 1 kHz, and show power spectra that cut off around 
500 Hz (Sunyaev & Revnivtsev 2000). 3) What is the mass? Without a binary companion, the mass of an accreting object is 
difficult to measure; however, this might be achieved by carrying out astrometry of background stars to look for gravitational 
distortion by the accreting object (Paczyhski 2001). 4) What do other parts of the spectra look like? Accreting black holes 
can produce relativistic radio jets (Mirabel & Rodriguez 1999), and photoionisation of the surrounding gas might result in 
observable infrared lines (Maloney, Colgan & Hollenbach 1997). 

There are several assumptions in our calculation that might affect the results. We have assumed that the accretion 
flow is one-dimensional, and that the sonic point is at the accretion radius. We have also assumed that the accretion is 
time steady; this may not be the case for higher accretion rates than we find for black holes accreting from the ISM (e.g. 
Grindlay 1978). If the black hole is moving slower than the sound speed at the accretion radius then the preheated region 
outside the accretion radius may have a chance to expand into the ISM, reducing the number density, and thus reducing 
the accretion rate. We have neglected mechanical feedback, which may occur if a jet or wind is developed. Recent work on 
non-radiating accretion flows indicates that strong winds can be formed which carry the bulk of the energy outwards as 
mechanical rather than radiative energy (Blandford & Begelman 1999, Igumenshchev, Abramowicz & Narayan 2000, Hawley, 
Balbus & Stone 2001). These authors argue that the accretion rate scales as M oc r. If the outer radius is taken as ta, 
then M(10r 3 ) ~ 10(i> / 'c) 2 M (r a) ~ 2 x W~ 7 v^g M{ta)- If the outer radius is taken to be the circularization radius, then 
M(10r 9 ) ~ 10~ 3 (M/9M Q )~ 2/3 [(i; 2 + c 2 ) 1/2 /40km s _1 ]~ 10/3 . Either of these circumstances will make black holes less visible, 
which in the formalism of this paper corresponds to a further reduction in the accretion efficiency. Magnetic fields will likely 
play an important role as they are amplified by flux-freezing, possibly heating the gas to virial temperatures (Igumenshchev 
& Narayan 2001) and can transport angular momentum via the magnetorotational instability once the circularization radius 
is reached, modifying the dynamics of the accretion flow. We have not explored the dependence of the efficiency on other 
parameters, which may result, for example, from changes in the gas density and angular momentum as a function of the 
accretion rate, so our extrapolation from estimated efficiencies of black-hole X-ray binaries may be too optimistic. The validity 
of these assumptions can be tested with 3-D radiation MHD simulations of a black hole accreting from an inhomogeneous 
medium, a daunting numerical problem. 

It is possible that a population of intermediate-mass black holes (IMBHs) exist with larger masses around 25OM0, the 
remnants of the first generation of star formation (Madau & Rees 2001). The larger masses of these objects would result in 
yet larger accretion rates and luminosities by a factor of ~ 10 3 , if they are distributed with same phase-space distribution as 
9Mq black holes. If we assume that ~ 10 6 IMBHs reside in our galaxy, then we find that the number of detectable objects at 
high fluxes may comparable to 9Mq black holes accreting at the same efficiency (Figure 3). However, if these objects reside 
in the halo or bulge of the galaxy, the number detectable will be decreased significantly. 
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